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Abstract. Tabe and Yokoyama found recently that the optical axis in a chiral monolayer of a ferronematic 
rotates when water evaporates from the bath: the chiral molecules act as propellers. When the axis is 
blocked at the lateral walls of the trough, the accumulated rotation inside creates huge splays and bends. 
We discuss the relaxation of these tensions, assuming that a single dust particle nucleates disclination 
pairs. For the simplest geometry, we then predict a long delay time followed by a non-periodic sequence of 
"bursts". These ideas are checked by numerical simulations. 

PACS. 61.30.Hn Surface phenomena: alignment, anchoring, anchoring transitions, surface-induced layer- 
ing, surface-induced ordering, wetting, prewetting transitions, and wetting transitions - 61.30.Jf Defects 
in liquid crystals 



1 Introduction 

The chiral compound R(OPOB) forms a bulk ferroelectric 



J w = br + PA\x 
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Fig. 1. Chemical structure of R(OPOB) 



liquid crystal. But it can also be deposited as a monolayer 
on water (or on water glycerol mixtures, which slow down 
the dynamics). Tabe and Yokoyama observed these lay- 
ers optically and found that the optical axis of this two- 
dimensional nematic was rotating continuously in the x-y 
plane . They interpreted this remarkable result in terms 
of an irreversible process: water evaporation from the sup- 
porting liquid is responsible for the molecular precession. 
They wrote down the basic Onsager relations required to 
describe the effect. They consist of 

a) Two fluxes: the rotation rate f2 and the mass flux 
of water evaporation J w . 

b) Two forces: the torque r (around a vertical axis) 
applied on the nematic molecules by external fields or by 
local distortions, and the chemical potential difference A[i 
between water in the liquid and water in the gas [2]. 

The entropy source per unit area of the monolayer is 
then 



In Eq. P|l. 71 is a classical friction coefficient for nematic 
rotations 3 and P is related to the permeability of the air 
layer above the Langmuir trough. The interesting coupling 
coefficient is b. 

Ref. £Q shows that in simple conditions (Z\/i is fixed 
and t ~ 0), fl is indeed proportional to Afi, and that the 
experimental values of b are of the right order of magni- 
tude (typically the rotation periods 2ir / ' fl are a few sec- 
onds). 

We are interested here in one particular feature - the 
angle 4> (labelling the direction of the optical axis in the 
horizontal plane) is forced by water evaporation to rotate 
in the monolayer. But at the lateral edges of the trough, 
the direction is expected to be anchored at the walls 0- 
This frustrated situation forces the internal clastic energy 
(proportional to |V</>| 2 ) to relax by successive "bursts". 
Our aim is to present a simple description of these bursts. 



2 Weak and strong anchoring 

The basic equation for the director angle (f>(x,y,t) in the 
monolayer can be derived from Eq. Q , incorporating elas- 
tic torques into r, and is 



dt 



7i 



n 



(4) 



TS = t[2 + AjiJ w 
and the Onsager relations are 

Q = 7 1 " 1 r + bA/j, 



(1) 



where Q is fixed because A[i is imposed. We have as- 
sumed that the splay and bend elastic constants are equal, 
Ki = K 2 — K. It will be convenient to define a diffusion 
coefficient 

7i 
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which we expect to be of order 10~ 6 -10~ 7 cm 2 /sec, de- 
pending on the water/glycerol fraction. 

Frustration originates from the anchoring at the lateral 
boundaries of the trough. We expect to have an energy per 
unit length at the boundary of the order — £/cos(0 — 0), 
where (3 define a preferred orientation (e.g. normal to the 
walls) When expressed per molecule, the anchoring en- 
ergy U is small compared to the nematic coupling (maybe 
10 times smaller), but it is very important. If we have 
a gradient V0| n normal to the wall, anchoring ruptures 
whenever the torque ifV0| n becomes larger than U . This 
corresponds to a critical value of the gradient 

*W = V<f>\ n = U/K ~ l/10a (6) 

where a ~ lnm is a molecular size. Thus we see that 
k s ~ 10 6 cm" 1 . Here we shall be mostly interested in the 
limit where the gradient |V0| is smaller than K SU rf, i-e., 
the strong anchoring limit. Then frustration effects in the 
interior of the monolayer are important. 



When this gradient reaches some critical value a 
"burst" may occur. But we have to be careful: this burst 
is a nucleation process for a new structure. Nucleation in 
nature is not controlled by the ideal (homogeneous nu- 
cleation) threshold. It is catalysed by a dust particle, a 
cosmic ray, or some other perturbation. A realistic for- 
mulation requires that we put a point-like "dust particle" 
somewhere in the trough. We assume that a burst will oc- 
cur when the magnitude of the gradient at this point |V0| 
reaches a critical value k. At all other points |V0| may 
possibly exceed n. 

Thus we start our discussion by looking at Fig. 3 and 
imposing that one dust particle is located at some point 
x = d and y — 0. We impose that d <C \J~D~i < L to be in 
the region described by Eq. Q. When the local gradient 
d(j)/dx located at the point x reaches the critical value k, 
a new scenario must start. 



4 The first burst 



3 The role of the dust particle 




Fig. 2. A nematic monolayer confined by two walls at x — 
and x = 2L, with boundary conditions <f> = at both ends. At 
t = 0, <f) = everywhere, dashed line. After time t has passed 
and under the influence of Eq. Q the profile evolves, solid line. 

We shall focus our attention on a simple geometry (see 
Fig. 2), where the monolayer is limited by two parallel 
walls separated by distance 2L. The boundary conditions 
are that — (modulus 2n) on each wall jS] . Assume, for 
instance, that at time t = we have = everywhere 
(Fig. 2, dashed line). At a later time t smaller than L 2 / D 
the angle <f> has rotated to <p = fit in a central region, but 
it drops to zero in a diffusion region of size (Dt) 1 / 2 near 
each wall (Fig. 2, solid line). In these regions we have a 
strong gradient 

50 fit n ft\ 1/2 

where the constant above, c = 2/y / 7r, is taken to be unity 
in subsequent calculations [Hj- 



Y 


A 




M 




R- 


dust particle 











o J 












) 




L X 


R- 


A; 









Fig. 3. Nematic monolayer confined between x = and x — L 
under the influence of driven rotations. The dust particle is at 
x — d and y = 0. When elastic stresses are too large, two point 
disclinations at A\ and A2 appear. 

The first breaking event occurs when d<f)/dx described 
by Eq. reaches the critical value k. This corresponds 
to a time t\ = n 2 D / Q 2 . We think of k values of order 10 5 
cm -1 (comparable to the inverse size of the dust particle). 
Thus t\ is of order 1 hour and the corresponding diffusion 
length is (Dti) 1 / 2 ~ 1 mm. Our solution after the burst 
implies a pair of disclinations A\ and Ai located at x = 
d and y — ±R(t), see Fig. 3. The initial value of R is 
comparable to the dust particle size. The <j) field is then 
described by 

= x ~~fr~~ + 0i = ^0 + 0i (8) 

where the first term corresponds to Eq. J7J) and 0i is due 
to the disclination, and is given simply by 

0i(Af) = a (9) 
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Fig. 4. (a) Ferro-nematic vector in the x-y plane from <f>\ of the 
pair disclination of Fig. 3. The two point defects are marked 
with filled circles, (b) orientation vector corresponding to the 
field <j> — <f>o + (j>\ given by Eqs. (|SJ and © . Empty and filled 
circles mark the initial and final location of the point defects 
in the simulation. Lengths are scaled by a = 10~ 7 cm. 



Here a is the angle A1MA2 defined as in Fig. 3, and 
M = (x,y) is the observation point, tfii satisfies the equa- 
tion V 2 (/)i = 0, as first noticed by F. C. Frank (see Fig. 
3). We assume that the distance 2R = A1A2 is small 
compared to the wall distance, R d and therefore the 
boundary conditions at the wall are satisfied without in- 
cluding images. Since V0i is opposed in direction to the 
unperturbed gradient V0Oj the disclination pairs reduces 
the elastic energy F c i 



F el = \k J (V^) 2 dxdy = E Q + E 1 + E 2 
where Eq is the unperturbed form, while 



£1 



can be reduced to a line integral along A1A2, 

E 1= XV^o^i] dl = -K\\7(j) \ 4ttR 

J A! 
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(12) 



where we introduced a cut along the (A\, A2) segment, 
and the discontinuity of 0i is [<f>i] — 2ir. Finally, E 2 is the 
energy of the disclination pair in an undistorted medium 



E 2 = ^K J {V fa) 2 dxdy = 2irK\n f— 



(13) 



Here a is proportional to the core radius of the disclina- 
tion. 

The force acting on A\ has an attractive component 
derived from Eq. (|13f) and a repulsive component derived 
from (O 

H = -\^i E i + 3») = {\Vcj> Q \ - (14) 

At the moment of burst |V(/>o| has its threshold value k. 
The regime of interest is R and the force is simply 

fl = 2ttkK. As is known for nematics [3], a disclination 
line moving in a constant force has a (nearly) constant 
velocity v 



R = h/C 



(15) 



where £ is a friction coefficient due to the rotations in- 
duced by the moving line 



(16) 



We conducted a numerical study where d<j>/dt = DV 2 4> 
was solved with different but convenient boundary condi- 
tions, namely <fi = at the left wall and <p — Ht at the 
right wall. The dust particle was assumed to be half-way 
between the walls. Figure 4 (b) shows the director in the 
x-y plane at some time after the first burst, i.e. after a 
disclination pair is created. It consist of a sum of the uni- 
formly changing field <po and the disclination given by <pi . 
Filled circles mark the defect cores while empty circles 
mark the cores just after the burst. 

Figure 5 (a) is a plot of <f> as a function of lateral coor- 
dinate x, for y = in Fig. 4 (b). The solid line corresponds 
to Fig. 4 (b) while the dashed line is just after the burst 
event. Clearly, the gradient <f>'{x) is lowered near the dust 
particle (x — 40). In Fig. 5 (b) we show the time evolution 
of the inter-pair distance R(t). In the simulation we have 
Rk 3> 1 and the nearly linear dependence of R on t is 
evident. 



(10) 5 The sequence of bursts 



We return now to the original problem, with an unper- 
turbed gradient described by Eq. J7J. It is easy to see 
that just after the first burst at t = t\ = k 2 D/Q 2 (and 
ignoring a factor of 1/c 2 ~ 0.8 in Eq. 10), the correction 
to the gradient near the dust particle is d(f)i/dx — —2/R, 
so that the gradient of (j) in the x-direction is reduced. 
After the burst, <f> increases with time. The second burst 
appears when the gradient becomes again equal to k, 



1 



(DtaJVa R t 2 
We introduce the notation 



A = (D/Q) 



1/2 



(17) 



(18) 
(19) 
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Fig. 5. (a) dashed line is a plot of <j> as a function of x just 
after the first burst. The dust particle is at x — 40 (a;-axis 
scaled by a = 10~ 7 cm). The gradient 4>'(x) is reduced near 
the particle. Solid line shows the profile at some later time, 
(b) plot of the defect location R(t) as a function of time. The 
near- linear behavior is in accord with Eqs. 1141 and 1151 . We 
used D = 10 -6 cm 2 /sec and k = 10 6 cm -1 . 



With this notation we find that 6\ — kX ~ 10 2 is a rather 
large number. Using R = Rt = Dnt and QX/Dk — 1/kX, 
Eq. I|17[) can be rewritten in the form 



= kX 



1 



Analysis of this equations reveals that 02 is rather close 
to B\ : treating n X as a large number one can obtain that 



kX + 



1 



V2kX 



(21) 



The second term on the right hand side is very small com- 
pared to the first one. What is the time difference be- 
tween the first and second bursts? If ft = 1 sec -1 then 



t% — t\ = 1.4 sec. Thus the bursts start to be frequent 
once we wait more than the long time t\ . 

Using the assumption k> 1 / R (meaning that the de- 
fects cf>i interact mainly with 4>$ but not with themselves), 
we get for the n'th burst event 



(Dtn) 1 / 2 £V*n-*»-i 1 



1 



— Ui-2 
1 



Or, written differently, 

8 n = kX + 



n— 1 

— y 

kx 4^ 



i 



9 2 - 9 2 



(22) 



(23) 



The solution for 62 obtained above [Eq. lt2*T|) ] suggests that 
the sequence of bursts has the form 6 n = kX + An a , with 
kX ^> An a . We put this form of 8 n and approximate the 
sum by an integral, 



An a 



1 

kX 



1 



dm 



2AkX n a - m° 



We change variable from m to p — m/n: 

,1-1/n ^ dp 
n/n 



An c 



2A(nX) 2 n a 



l-p° 



(24) 



(25) 



This equation is satisfied if a = 1/2. In this case the in- 
tegral on the right hand side is equal to | ln(n — 1) + 
~ ln((2n + 1) /(n + 1)), and can be approximated by unity 
because of the logarithmic dependence on n. In this ap- 
proximation we find that A = ^gjjj. In summary, the 
burst sequence is given by 



9 n = nX+^^-n 1 ' 2 
V2kX 



(26) 



The time difference between two successive bursts increases 
as n 1 / 2 , 



1 hi 

and is of the order of few seconds. 



(27) 



( 20 ) 6 Conclusion 



We have studied a ferro-nematic monolayer under the in- 
fluence of water evaporation. This evaporation introduces 
an external torque tending to rotate the individual molecules. 
Frustration should occur because the molecules at the 
edges of the Langmuir trough are anchored to the walls. 

For a simple geometry, we show that a dust particle or 
some other disturbance can initiate a pair of disclinations, 
thus relieving the stress near the particle. The distance 
between the two disclination points increases linearly with 
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time, until the local strains near the particle again reach 
their threshold. At this time two more disclinations are 
born, also moving away from each other, and the process 
continues. 

What would happen in a more realistic geometry, for 
instance with a circular trough, and anchoring condition 
which imposes a director normal to the walls? In this case, 
the ideal starting configuration at t = would have a 
disclination point at the center of the trough. When we 
switch on the evaporation process at t — 0, the align- 
ment around this point would rotate uniformly in time: 
4>{x, y, t) = (j)(x, y, 0) + fit. But near the walls, in a zone 
of width {Dt) 1 / 2 , we would find strong gradients, and dust 
particles could initiate certain "ladders of bursts" . 

There are few practical complications to be expected: 

a) The initial state will usually involve many disclination 
points rather than one. We still expect a uniform rota- 
tion of the pattern in the central region - this probably 
corresponds to the experiment of Ref. pQ. One could re- 
move the disclination from the central region by preparing 
the sample under a horizonal magnetic field H , and then 
switching off the field at t = 0. 

b) Our description of the nemato-hydrodynamics was sim- 
plistic jnj. With a fuller description we might expect cer- 
tain backflow effects - the dust particles might move in 
the x-y plane. 

c) The width (Dti) 1 / 2 of the diffusion layer (at the on- 
set of the ladders) is expected to be rather small, about a 
millimeter. In this region, the monolayer need not be flat, 
as a meniscus will be created near the walls. 

A final remark: there is a superficial analogy between 
rotations due to evaporation, which we discuss here, and 
the rotations induced in a nematic by a rotating magnetic 
field 3 . But there is a deep difference. When the field is 
present, the rotational symmetry is broken, and the de- 
fects generated by frustration are domain walls. But in 
our case, rotational symmetry is preserved and the rele- 
vant defcts are disclination lines. 



7 Appendix- burst sequence in one dimension 

We now discuss the series of burst events in the one- 
dimensional case, namely <j) depends on one variable only, 
x. The time t\ where the first event occurs is given by 



fit 1 



(Dh) 1 / 2 



(28) 



In our prescription, the defect is created in a "delta func- 
tion" manner, i.e. one molecule changes its orientation. Af- 
ter some time, this single molecule affects its neighborhood 
in a diffusive manner. Hence we write that the slope is re- 
duced by a value ^Tr/D{t — ti) exp ^— (x — d) 2 / \j2D{t — t x 

where d is the particle location. At the second burst event, 
the total slope (j)'{d) is again equal to the threshold value 
k, and hence we have at t 2 



Similarly, the n'th event is given by 

fltn 1 7T / 1 



(Dt n y/ 2 




y/^ri t-n— 1 , 



(30) 



Using 6 n = (fit n ) 1/2 and A = (D/fi) 1 / 2 , we rewrite the 
last equation as 



2 _ 02 



(31) 



We continue to solve this relation assuming that 9 n = 
An a + Ak, and that on the right hand side of Eq. (|31|l . An 
is larger than the sum. We also approximate the sum by 
an integral, and find that 



An a = 



rl/2 



\/2\KAn a 



dm 



1- (? 



(32) 



This equation is satisfied when a = 2/3: in this case the 
variable in the integral can be changed from m to p = nm, 
and the integral becomes equal to 37r/4. In summary, 



9 n = An a + Xk 
n = 2/3 



(33) 



A = 



1 



-7T 



3/2 



fiU 



(Dt 2 y/ 2 ^D(t 2 -h) 



(29) 



The above sequence of events could be relevant for a 
one dimensional superconductor rod. If the rod is thin, in 
each cross-section the phase is constant. When a current 
is driven into the system, this phase changes continuously 
along the rod. When the phase gradient is too large, a 
similar defect should be created, nucleating most probably 
around a mechanical "weak spot" . The detailed nucleation 
process was discussed long ago by Langer and Ambegaokar 
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